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ABSTRACT 

Using observations obtained with the LOw Fequency ARray (LOFAR), the Westerbork Synthesis Radio Telescope (WSRT) and 
archival Very Large Array (VLA) data, we have traced the radio emission to large scales in the complex source 4C 35.06 located in 
the core of the galaxy cluster Abell 407. At higher spatial resolution (~ 4"), the source was known to have two inner radio lobes 
spanning 31 kpc and a diffuse, low-brightness extension running parallel to them, offset by about 11 kpc (in projection). 

At 62 MHz, we detect the radio emission of this structure extending out to 210 kpc. At 1.4 GHz and intermediate spatial resolution 
(~ 30"), the structure appears to have a helical morphology. 

We have derived the characteristics of the radio spectral index across the source. We show that the source morphology is most likely 
the result of at least two episodes of AGN activity separated by a dormant period of around 35 Myr. The outermost regions of radio 
emission have a steep spectral index (a < -1), indicative of old plasma. We connect the spectral index properties of the resolved 
source structure with the integrated flux density spectral index of 4C 35.06 and suggest an explanation for its unusual integrated 
flux density spectral shape (a moderately steep power law with no discernible spectral break), possibly providing a proxy for future 
studies of more distant radio sources through inferring their detailed spectral index properties and activity history from their integrated 
spectral indices. 

The AGN is hosted by one of the galaxies located in the cluster core of Abell 407. We propose that it is intermittently active as it 
moves in the dense environment in the cluster core. In this scenario, the AGN turned on sometime in the past, and has produced the 
helical pattern of emission, possibly a sign of jet precession/merger during that episode of activity. Using LOFAR, we can trace the 
relic plasma from that episode of activity out to greater distances from the core than ever before. 

Using the the WSRT, we detect HI in absorption against the center of the radio source. The absorption profile is relatively broad 
(FWHM of 288 kms“^), similar to what is found in other clusters. The derived column density is Nm ~ 4 • 10^*^ cm“^ for a 
Tspin = 100 K. This detection supports the connection - already suggested for other restarted radio sources - between the presence 
of cold gas and restarting activity. The cold gas appears to be dominated by a blue-shifted component although the broad HI profile 
could also include gas with different kinematics. 

Understanding the duty cycle of the radio emission as well as the triggering mechanism for starting (or restarting) the radio-loud 
activity can provide important constraints to quantify the impact of AGN feedback on galaxy evolution. The study of these mechanisms 
at low frequencies using morphological and spectral information promises to bring new important insights in this held. 

Key words, galaxies: active - radio continuum: galaxies - radio lines: HI - galaxies: individual: 4C 35.06 


1. Introduction 

For some time now it has been recognized that the active phase 
of radio-loud active gala ctic nuclei (AGN) can be recurrent (see 
ISaikia & Jamrozvll2009L for a review). Therefore, this activity 
can have a repeated impact on the AGN host galaxy and its envi¬ 
ronment, starting from the initial phase of the radio source’s life 


dO’ Deall2002h on to its later stages (iMcNamara et al.ll200Ql) . The 
morphol ogy of a (recurrent) mdio source can vary sign ificantly 
(see e.g. ISchoenmakers et al.ll200(il: iMurgia et al.lf201 lb . There¬ 
fore, only a combination of radio morphology and radio spectral 
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inde43 analyses can reliably tell us whether a radio source had re¬ 
curring activity, thereby affecting its host galaxy more than once. 
A typical radio-loud AGN in its active phase (which can last up 
to several 10^ years) is characterized by the presence of a radio 
core and extended lobes with or without visible jets. When the 
AGN shuts down, emission from the core disappears for the du¬ 
ration of the "off" phase. The large scale radio lobes might still 
be replenished with particles if the information of the core shut¬ 
down has not had enough time to propagate out. However, once 
the supply of energetic particles coming from the core stops, the 
lobe plasma starts to age producing a so-called relic structure and 
the radio source enters a fading phase in its life-cycle. It is possi¬ 
ble that during this time the AGN re-activates, in which case we 
can again observe an active radio core (possibly with new jets / 
lobes expanding outward into the ISM / IGM) co-existing with a 
fading radio relic on larger scales. 

The energy loss mechanism for relic radio plasma is mostly 
through synchrotron radiation and inverse Compton (IC) scat¬ 
tering off of cosmic microwave (CMB) photons (with possible 
adiabatic expansion related losses). Its observational signature is 
a steepening of the radio spectrum of the emission regions (syn¬ 
chrotron radiation and IC losses deplete high-energy particles 
faster than low-e nergy ones), maki ng them relatively brighter at 
low frequencies (iKardashevI 1 1 962h . Thus, the relic radio emis¬ 
sion is typically faint and detecting it requires low-frequency 
observations with high sensitivity. 

This work is part of a larger study where we investigate 
radio sources with restarted activity expanding our studies to 
low frequencies, using the LOw Frequency ARray (LOFAR, 
IVan Haarlem et al.ll2013l) . Our aim is to describe in more detail 
their activity history (duty cycle). We also benefit from the ad¬ 
vantage of having sufficient spatial resolution to be able to study 
resolved source structure at low frequencies. We relate our re¬ 
sults to the integrated flux density spectral properties of these 
sources in an attempt to extrapolate our results to the high red- 
shift source population (in which case we would have poorer 
spatial resolution). 

In this work, we investigat e the complex radi o morphology 
of 4C 35.06 (B2 0258+35B: iBondi et alJ [To^ . iLiuzzo et all 
(I 2 OIOI) did a study of UGC 2489 using the Very Long Baseline 
Array (VLBA) and they detected an extended source on a milli 
arcsec scale with a core flux density of ~ 2.6 mJy, stating that 
the VLBA detection coincides with one of the cores of the cD 
galaxy UGC 2489 (Figure IB. UGC 2489 represents th e core of 
the Abell 407 galaxy cluster (ISchneider & Gunnil 19821 see Fig- 
ure[T]). The radio luminosity of 4C 35.06 is L\ agy\j = 2.5 • 10^"^ 
WHz~ ^. This value is consistent according to iBest & HeckmanI 
(l2012h with being a low excitation radio galaxy (LERG) as 
confir med by the low-reylution optical spectm of the host 
galaxy (ISchneider & Gunnil 1 9821: ICrawford et ^1 19991) showing 
no emission lines. 

In the VLA maps we can see two inner radio lobes spanning 
31 kpc and a diffuse, low-brightness extension running parallel 
to them, oriented southeast to northwest (see Figure[T]). The puz¬ 
zling structure hints at the possibility of multiple AGN activity 
episodes combined with motion of the host core within the cD 
galaxy, making this source particularly interesting for studying 
AGN recurrent activity and the triggering of the active phases. 

We made use of images available in the Jansky Very Large 
Array (VLA) archive. Furthermore, we took advantage of the 


^ We use the 5* oc y® definition for the spectral index a throughout this 
writing. 
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Fig. 1. Optical image of UGC 248 9 taken from the Sloan Digital Sky 
Survey fSDSS. lAihara et aOLOllh showing its multi-component core 
structure. The red intensity overlay shows a 4885 MHz VLA archival 
image. In yellow we over-plot contours from a 1415 MHz VLA archival 
image. Contours are log scaled, there are 10 contour levels going from 
4cr to 500cr. Details about the VLA images can be found in Table[T] The 
stellar halo around the cores is roughly 40 kpc in extent. The position 
of the core hosting the AGN is shown with a black cross. 


fact that this source was included in one of the fields observed 
by the LOFAR during commissioning tests. 

In addition, we also present new Westerbork Synthesis Ra¬ 
dio Telescope (WSRT) observations obtained to investigate the 
presence of cold gas (HI) in this object. As mentioned above, 
one of the open questions is the relation between the gas and the 
fueling of multiple phases of activity. Atomic neutral hydrogen 
can, in principle, provide such fuel. 

The paper is organized in the following way. Section [2] dis¬ 
cusses the data reduction procedure, while in Section [3] we 
present the LOFAR image and quantities derived from it. We 
discuss the implications of our results, in combination with re¬ 
sults of previous studies in Section |4l and we conclude with a 
summary in Section [5l 

Throughout this paper we assume a fiat A cold dark matter 
(ACDM) cosmology with Hq = 73 km s“^ Mpc“\ Qa = 0.73 
and Qm = 0.27. At the distance of 4C 35.06, using z = 0.046276 
(corrected to the referen ce frame defin ed by the CMB), this re- 
suits in 1 " = 0.872 kpc (IWrighlll2006l) . 


Table 1. Summary of the images of 4C 35.06 used in this work. 


Instrument 

y [MHz] 

Beam size 

Noise [mJy/b] 

Ref. 

LOFAR 

61 

51"x43" 

35 

1 

WSRT 

1360 

30"x 13" 

0.18 

1 

VLA (B) 

1415 

4'.'8 x 4'.'1 

0.29 

2 

VLA (C) 

4885 

378 

0.04 

3 


Referenc es. (1) This work; (2) jGregorini et al.ll 19881 . VLA archive im¬ 
age); (3) dBondi et al.ll 19931 . VLA archive image). 
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2. Data available, new observations and data 
reduction 

2.1. LOFAR 

The observations using LOFAR were carried out in the night of 
October 8, 2011 using the low band antennas (LBA) at 62 MHz 
centered on B2 0258+35A observed as part of the LOFAR com¬ 
missioning efforts. 4C 35.06 is located about half a degree to the 
north from the phase center. Table [2] summarizes the observa¬ 
tional configuration. 

Table 2. LOFAR observational configuration 


No. of sub-bands 

28 

Central Frequency [MHz] 

62 

Bandwidth [MHz] 

5.6 

Integration time 

1 second 

Observation duration 

6 hours 

Polarization 

Full Stokes 

No. of LOFAR stations used 

32 

UV range 

20+ - 4000+ 


The data were flagged to remove radi o frequency interfer- 
ence (RFI) using the AO flagger package (lOffringa et al .11201^ 
and averaged to a frequency resolution of 192 kHz per chan¬ 
nel and a time resolution of 10 seconds per sample. Calibra¬ 
tion was performed using the Black Board Selfcal (BBS) pack¬ 
age using the complete bandwidth. The set up of the observa¬ 
tions did not include observations of a flux calibrator, therefore 
the initial calibration model was deri ved from a VLSSR image 
foll owing the procedure desc ribed in IVan Weeren et alT (l2012b 
and iDe Gasperin et~^ (l2012h . The image from which the initial 
model was derived had a size of 8° x 8° to encompass the full 
width at half maximum (FWHM) of the LOFAR station (pri¬ 
mary) beam, and to ensure that any brighter sources beyond the 
FWHM are included in the model. The model components were 
specified as having a spectral index of a = -0.8. During the 
calibration, we have solved for the complex valued station gains 
and for the differential total electron content (dTEC) over the full 
bandwidth. 

Imaging was done using the Common Astronomy Soft- 
ware Applications pa ckage (CASA: iMcMullin et al.ll2007l) . us¬ 
ing multi-scale clean (ICornwellll2Q08l) . We imaged bv selecting 
the data 'mUV distance below 20 km, limiting ourselves to sta¬ 
tions having satisfactory calibration solutions, avoiding as much 
as possible calibration errors induced by the ionosphere, while 
retaining enough spatial resolution. We weighted the baselines 
used for imaging using Briggs weights with the robust param¬ 
eter set to -0.1 (lBriggs|[T995l) . We have used multi-frequency 
synthesis (assuming no frequency dependence of the sky model) 
imaging over the entire bandwidth. 

After obtaining the initial image, we have performed a cor¬ 
rection for the primary (station) beam shape to ensure a correct 
flux scale over the FoV. The correction was applied by deriving 
the average LOFAR station beam, using CASA to image (with 
the same settings as were applied in the imaging) a grid of BBS 
simulated point sources covering the LOFAR FoV. The recov¬ 
ered flux densities were used to derive the primary beam which 
we use to correct the image. 

Before the final imaging, three self-calibration runs were per¬ 
formed by iterating over the procedure described above. We have 


^ VLSS dCohen et al.l 120071) is the VLA Low frequency Sky Survey 
carried out at a frequency of 74 MHz 


reached an image noise of about 35 mJy beam“^ measured away 
from bright sources. The size of the restoring beam is 5 l"x 43". 
We have checked the flux scale across the final image, by com¬ 
paring the flux densities of sources in our LOFAR image with 
their values in existing surveys. Compact sou rces with S IN > 10 
were extracted using the PyBDSM package (iRamanujamll2Q07l) 
and their flux density at 62 MHz compared with the extrapola- 
tion from the VLSS, WENSS0 and NVSS0 catalogs. The derived 
LOFAR flux densities were estimated to be accurate to ~ 10% 
level. 

A check of the positional accuracy of the LOFAR image was 
performed by plotting the differences between positions of the 
sources extracted from the LOFAR image and matched sources 
from the NVSS survey. Again, we used PyBDSM to extract 
point sources from the LOFAR image. Positional differences for 
source matches between the WENSS and NVSS catalogs were 
used as control. On average, the derived position offset for point 
sources in the vicinity of 4C 35.06 extracted from the LOFAR 
image and the corresponding NVSS catalog sources is less than 
5" without any obvious systematics. Thus, for the purpose of 
this study, we concluded that there is no position offset in the 
LOFAR image. 

2.2. HI - 21cm and continuum from the WSRT 

WSRT observations to search for HI in 4C 35.06 were per¬ 
formed at three epochs (2013 April 8 and 11, June 12) for 8, 
7 and 6 hours respectively. The observing band was centered 
on 1356.8 MHz, the frequency of HI for z = 0.0464, the red- 
shift of 4C 35.06. The total bandwidth was 20 MHz, corre¬ 
sponding to a velocity range of ~ 4422 kms“^ and was cov¬ 
ered with 1024 channels (dual polarization). In order to increase 
the signal-to-noise ratio, four channels were binned resulting 
in a spectral resolution of 17 kms“^ in the output data cube. 
The data reduction was perform ed using the MIRIAD package 
(ISault. Teuben & Wright! [19^ . After flagging, bandpass and 
phase calibration were performed. The continuum emission was 
subtracted, using the task uvlin, by making a second-order poly¬ 
nomial fit to the line-free channels of each visibility record. 

Cubes with different weighting were obtained. For an inter¬ 
mediate weighting (Briggs with the robust parameter set to 0.4), 
the final r.m.s. noise was 0.2 mJy beam“^chan“^ for a velocity 
resolution of 15 kms“^ and spatial resolution of 34" x 14". In 
addition to the line cube, a continuum image was obtained using 
the line-free channels. The peak of the continuum emission was 
221 mJy and the r.m.s. noise was 0.18 mJy beam“^. 


3. Results 

3.1. Radio continuum structure from WSRT and LOFAR 

Our LOFAR and WSRT images of 4C 35.06, as well as data 
derived from the archival VLA images are presented in Fig. |2l 
Looking at the morphology traced by the different images at dif¬ 
ferent resolution, we can note several things. 

Firstly, Figure ^indicates that the emission parallel to the in¬ 
ner radio lobes is connected to the more distant regions through 
a helical structure (best seen in the left panel of Figure |2] as a 
bent tail of emission extending to the northwest from the cen¬ 
tral region, and then turning to the south; we note that there 

^ WENSS stands for the Westerbork Northern Sk y Survey carried out 
at a frequency of 325 MHz jRengelink et al.ll 19971) 

NVSS stands for the NRAO VLA S ky Survey carried out at a fre¬ 
quency of 1.4 GHz jCondon et al.ll 19981) 
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Fig. 2. Left: WSRT 1360 MHz continuum image of 4C 35.06. Overlaid in red are contours taken from the same image, with levels (-5,15,150)cr, 
while plotted in black are 4885 MHz VLA contours of the central region with contour levels of (-5,50,100,200,350)cr. The WSRT beam is 
shown in the lower left corner. Middle: The 4C 35.06 radio source as seen by LOFAR at 62 MHz, represented by red contours (33 levels spanning 
(-3,183)cr) over an inverted grayscale SDSS image of the central region of the Abell 407 cluster. The peak flux density in this map is 9.1 Jy. 
Regions in the LOFAR map are labeled A through E. Overlaid in blue are the VLA 4885 MHz image contours. Right: LOFAR 62 MHz (red) and 
VLA 1415 MHz (blue) contours overlaid over an inverted grayscale SDSS image of the cluster core (UGC 2489). The location of the core hosting 
the AGN is indicated with a white star. Arrows show the extent and the location of the regions of interest in different images. The relevant beam 
sizes and noise levels of the images are given in Table [T] 


is also a similar almost symmetric extension to the southeast). 
The WSRT image has sufficiently high resolution to enable us 
to justify the claim that the helicity is real and not an artifact of 
the elo ngated beam. It has recently been imaged by iBiiu et al] 
(12014l) using the Giant Metre Wave Radio Telescope (GMRT). 
Observations using LOFAR at a frequency of 150 MHz and 
higher resolution than described in this work also show the heli¬ 
cal substructure (Marisa Brienza, private communication). The 
helicity might indicate precession in the radio jet which has 
deposited the plasma, triggered b y accretion disk instabilities 
dKunert-Bajraszewska et al. I [20101) or maybe a possible sign of 
a past black hole merger of two cores in UGC 2489. 

The zoom-in of the central region of 4C 35.06 in the LO¬ 
FAR image (labeled C) is shown in the right panel of Figure [2] 
along with the same region as observed by the VLA at higher 
frequencies. 

The source region labeled E is noticeable as a barely visible 
surface brightness enhancement in a VLSS ima ge of 4C 35.06, 
and it i s detected at a 2cr level in a VLA image bv lGregorini et al.l 
(Il988l) . It is now unambiguously detected in the LOFAR image. 
Its projected distance from the central region (C) is ~ 210 kpc. 

The region labeled A was previously detected by the VLSS 
(blended with the rest of the source) and partially by the 
WENSS, however our LOFAR map detects it with better reso¬ 
lution and allows us to clearly morphologically separate it from 
other regions of 4C 35.06. 

The second interesting result is that, as shown in the right 
panel of Figure [21 the peak flux of 4C 35.06 as seen by LO¬ 
FAR is offset by 10" to the northwest with respect to the radio 
lobes seen in the VLA maps. This is a sufficiently large value to 
be considered a real effect and not a systematic offset (see Sect. 
12.11) . We interpret this difference between the peak flux density 
position and overall source position as intrinsic in the case of the 
extended source 4C 35.06 and indicative of its (spectral) proper¬ 
ties. 


3.2. HI gas 

Atomic neutral hydrogen is detected against the central region 
of 4C 35.06. Figure [3] shows the HI profile (obtained from the 
combination of all WSRT data) extracted against the brighter 
(western) lobe of the central source (Figure [T]). The profile is 
obtained after Hanning smoothing and averaged over a region 
encompassing the western inner lobe. Given the limited spatial 
resolution of the WSRT HI observations, we can only say that 
HI is seen in from of the western lobe which, in fact, includes 
the radio core of the source. 

The profile is complex, with a broader component with 
FWHM of 288 kms“^ centered at 13858 kms“^ The peak ab¬ 
sorption is ~ 0.8 mJy which corresponds to a low optical depth 
ofr = 0.0036. From the parameters of the absorption we can es¬ 
timate a column density of A^hi ~ 4 • 10^^ cm“^ for a Tspin = 100 
K. 

The HI profile is broader than that typicall y found in ra¬ 
dio galaxies (iGereK Morganti & Qosterlool2014l) . but resembles 
what is observed in some other clusters; we discuss this in more 
details in Section [4Hj While at uniform weighting the profile is 
dominated by two relatively narrow components, the broad com¬ 
ponent becomes more prominent at lower spatial resolution, sug¬ 
gesting that it belongs to a diffuse distribution of gas. However, 
given the low S/N of the detection this result will need to be ver¬ 
ified with deeper observations. Compared to the systemic veloc¬ 
ity of the host galaxy (14120 km s“^), the peak of the HI profile 
appears blueshifted (see Figure [3]) . However, the total width of 
the profile encompasses the range in velocity of the galaxies in 
the cluster core. 

In addition to the absorption profile, there are at least three 
detections of HI emission in the WSRT FoV. All three have 
an optical counterpart in the SDSS. One of these detections 
is identified with UGC 2493, a (disturbed, possibly merging) 
galaxy pair having a redshift of z =Q.0437, which m akes it a 
member of Abell 407 according to IChen et al.l (Il998l) . It is lo¬ 
cated at a radial distance of 210 kpc from the Abell 407 cluster 
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Fig. 3. HI absorption profile in 4C 35.06 (Hanning smoothing applied). 
The full vertical line indicates the systemic velocity of 4C 35.06 host 
galaxy, while the two vertical dashed lines show the highest and lowest 
velocity of the galaxies in the cluster core. 


core (UGC 2489). The other two detections are located at co¬ 
ordinates: J030219+354601 (HI disk associated with a edge on 
galaxy) and J030217-f355533. 

3.3. Spectral analysis 

Following the work of iMurgia et al.l (1201 ll) . we investigate the 
characteristics of the integrated spectrum of the source and then 
we turn to the spatially resolved analysis of the spectral index 
with the aim of getting a detailed handle on the past and present 
activity history. 

We have compiled integrated flux density measurements of 
4C 35.06, as shown in Table[3] The flux density values are taken 
from the literature, or measured from archival images, with com¬ 
parable sensitivity to the extended emission. They are plotted in 
Figured] along with our LOFAR flux value. 

The best power-law flt to the integrated flux density data 
points has a spectral index of a = -1.01 + 0.02. Although this 
is not an extremely steep spectral index, it is steeper than what 
is usually seen in the lobes of radio galaxies (a ~ -0.7 + 0.1). 
The integrated spectrum does not show a spectral break. 

This indicates that we average over regions of different ages 
having different intrinsic brightness. The sum of those spectra, 
each having spectral breaks at different frequencies, produces 
the remarkably straight integrated spectrum (see the discussion 
below for details). 

High-redshift radio sources that will be detected in future ra¬ 
dio surveys will be (in most cases) unresolved. Therefore, it is 
important to benchmark the information obtained from the inte¬ 
grated flux density spectral studies against the one for the spa¬ 
tially resolved studies. 

To expand on the different phases of radio-loud activity in 
this object, we used our LOFAR and WSRT observations, as 
well as archival VLA images to distinguish between the spec¬ 
tral index of the inner regions of 4C 35.06, and the outer regions 
of extended radio emission. The resulting spectral index images 
are shown in Figure [5| 

The maximum UV range of the LOFAR data set is 4kA, while 
that of the WSRT is 7 kT x 13 kA. The UV coverage within the 
range specified for the WSRT is uniform, while the LOFAR data 
set covers all of the spacings even though it has gaps due to the 


Table 3. Data used in Figure]?] 


Source v [MHz] S int [Jy] S c [Jy] 


Clark Lake^ 

26.3 

LOFAR^ 

61.6 

VLSS" 

74.0 

Cambridge^ 

81.5 

4C^ 

178.0 

MIYUN^ 

232.0 

WENSS^ 

327.0 

B2^ 

408.0 

WSRL 

1360.0 

NVSS^’ 

1400.0 

VLA^ 

1415.0 

VLA^ 

1490.0 

BKB77"^ 

2700.0 

GB6^ 

4850.0 

VLA^ 

4885.0 


47.0 + 6.0 

17.04 + 3.41 9.1 + 1.9 

13.92 + 0.34 

13.80+1.8 

5.01 + 0.87 

4.43 + 0.05 

3.27 + 0.01 1.6+ 0.4 

2.30 + 0.20 
0.76 + 0.03 
0.750 + 0.002 

0 . 6 + 0.1 

0.710 + 0.003 
0.37 + 0.02 
0.23 + 0.03 

0.20 + 0.01 


Notes. IVerkhodanov et al.l (120091) : Iviner & EricksonI (Il975h This 
work ICohen et ^ ( 20071). CLEAN bias correction a pplied 
(+ 0.7 Jy) 

IPilVino 


(d) 


menetal.! (izuu/l). CLJiiA J N bias correction a pplied 
IVerkhodanov et al.l (120091) : iReadhead & Hewish] il974l) 


to 


X 1.067 


(corrected according 

[ to iHelemboldt et al.l 

2008 

) jVerkhodanov et al. 

2009 

); IZhans et al.l 

(119971) iRenselink et al. 

(119971) <« iCollaetal. 

1973 

), flux set to 

lBaarsetal.1 (Il977l) x 1.091 (corrected accord- 


VLA B array. Table [T] V L A C array. IVerkhod anov et af 
(m) IVerkh odanov e t~^ d^2009l): iBridl e. Kestev en & Brandid 
IVerkhodanov et al.l (120091) : lGregorv & CondonI d 19911) VLA C ar¬ 

ray, Tabled] 


( 1998 ) 

(2009) 

]V7% 


observation length. We have smoothed the WSRT image to high¬ 
light the shorter spacing’s. The VLA archival images were ob¬ 
tained from data sets with matching UV ranges and uniform UV 
coverage. 

We see that the inner region where the core and the two 
brighter lobes are located shows a typical spectral index in the 
range of cr ~ -0.6 to cr ~ -0.8, indicative, as shown by the 
morphology of the inner double-lobed source, of an active re¬ 
gion. However, there is a spectral steepening (a ~ -1.2) in the 
extended emission just to the northwest of the inner lobes (Figure 
[5] left panel). This is the region of the source where we detect the 
peak flux density in our LOFAR map (see Figure [2] right panel), 
in agreement with LOFAR being more sensitive to regions hav¬ 
ing steeper spectrum. Thus, this region could indeed represent a 
region of older radio plasma. 

We used lower resolution images to derive a spectral index 
map for the larger scale regions of the source. The steepening 
seen at higher resolution in the structure running parallel to the 
inner brighter lobes, is also visible and continues out to larger 
scales. This can be seen in the right panel of Figure [5] Going to 
the northwest, the spectral index remains at a value of around 
-1.2 up to ~ 50 kpc from the central region, while on the south¬ 
eastern side it steepens already to about -1.4 quite close to the 
central region. This trend suggests that the outer, extended, re¬ 
gions of the source are composed of aged plasma. 

We have fitted a synchrotron aging model including a period 
of activ ity and a period during which the AGN was switched ofl T 
(KGJP; iKomissarov & Gubanovll 19941: IVan Weeren et al.ll2009l) 
to the integrated flux density measurements at four different fre¬ 
quencies (images listed in Table [B for the central region of the 
source (labeled C in Figure [6]). We have ch osen the KGJP model 
(Komissarov-Gubanov prescription with a lJaffe & Perolal[T97^ 
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Fig. 5. Left: spectral index map of the inner lobes of 4C 35.06 using the archival VLA images listed in Table[T] Overlaid are VLA 4885 MHz 

contours (-5,18,27,36,81)cr. Right: spectral index map of the extended emission of 4C 35.06 using the LOFAR and WSRT images from 

this work. Overlaid are 33 contour levels from the LOFAR map spanning the range (-3,183)cr. The beam sizes are given in the lower left corners. 



Fig. 4. Integrated flux density values for 4C 35.06 based on the data 
listed in Table [3] (circles) and its central region (stars). The LOFAR data 
points are in red. The green line represents a power law fit, while the 
blue dashed line gives the synchrotron aging model fit for region C. 


loss term) since region C contains the radio emission of the cur¬ 
rently active AGN as well as an older emission as is evident from 
the spectral index mapping. We have used an injection spectral 
index of ao = -0.8. The equi-partition value for the magneti c 
field was calculated to be 5.34/iG (according to lMilevlll980l) . 
Our best fit model gives a source age of 4 = 4n + 4fF = 1.5 + 10 
Myr. Here, 4 signifies the total source age, comprised of the time 
the source was active (4n) and the time elapsed since the AGN 
shutdown (4ff). The model supports the observation that the cen¬ 
tral region is the youngest region of the source. We show the fits 
in Figure m 


We have also determined the spectral index of the outermost 
source region labeled E, as well as that of region labeled A by 
measuring the flux density of these regions as seen in the LOFAR 
image. We have adopted the measurement procedure outlined by 
iKlein & EmersonI (1 19811) . The measurements were done by inte¬ 
grating the flux density over measurement areas in the regions 
A and E where the surface brightness was above the 3 cr level in 
the LOFAR, VLSS, WENSS and NVSS images, as well as in our 
continuum WSRT image and thus calculating the spectral index. 
The derived values are given in Tabled 

We plot the data from Table |4| in Figure [6] and lit a power- 
law to the data points for region A. The spectrum of region E is 
poorly constrained, but we note that it is ultra steep (a < -2). 
Region A has a flatter spectrum (a = -1.36 ± 0.09) in line with 
what we expect from the spectral trend of the inner regions. More 
sensitive observations at higher frequencies (with sufficient res¬ 
olution) are needed to determine the exact spectral index of the 
outermost regions. In Figure[6]we fit a second-order polynomial 
to the measurements to illustrate the aged spectra. 

4. Discussion 

Figure [T] shows the radio emission around the galaxies compris¬ 
ing the Abell 407 core (cataloged as UGC 2489). We can see that 
there are two lobes detected by the VLA, and diffuse emission 
that extends to the northeast. Going farther out, this emission 
extends into a helical extension as seen in the grayscale WSRT 
image shown in the left panel of Figured 

The extended emission to the northwest continues out to an 
even larger distance than previously known, i.e. ~ 210 kpc from 
the central region of 4C 35.06. It is labeled as region E in the 
LOFAR image (see Figure O. 

The overall source morphology is suggestive of (at least) two 
distinct episodes of radio activity in 4C 35.06. The extended 
emission is a remnant of the previous active phase, whilst the 
twin lobes detected by the VLA show the current activity cycle. 

Using the morphological and spectral results described 
above, we propose that 4C 35.06 can be interpreted as a source in 
which we are seeing (at least) two subsequent episodes of AGN 
activity. 
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Fig. 6. Left: The dashed blue line represents a second-order polynomial fit for region A. Middle: The 4C 35.06 radio source as seen by LOFAR, 
represented by red contours over an inverted grayscale SDSS image of the central region of the Abell 407 cluster. Overlaid in blue are the VLA 
4885 MHz image contours. Black circles show the measurement areas in regions A, C, and E. Right: Second-order polynomial fit (dashed blue 
line) for region E. The data points which are upper limits are represented with empty markers and downward pointing arrows. They are not taken 
into account in the fitting. 


Table 4. Measured flux density values for the regions labeled A and E 
in the LOEAR image. 


y [MHz] 

[Jy] 



Region A 

a 

61.6 

0.77 + 0.16 

-1.36 ±0.09 

74.0 

0.51+0.12 


325.0 

0.075 + 0.005 


1360.5 

< 0.0023 



Region E 

^62 

61.6 

0.96 + 0.20 

< -2^ 

74.0 

0.32 + 0.12 


325.0 

< 0.022 


1360.5 

< 0.0015 



Notes. Region E is detected in the LOEAR and VLSS images, and 
the other measurements are only upper limits. Hence, we only show an 
illustrative polynomial fit in Eigure 0 but do not attempt to derive a 
spectral index. 


To explain the complex radio morphology as well as the 
spectral indices we have derived, we propose that the AGN is 
intermittently active as it moves in the dense environment in the 
cluster core. In this scenario, the AGN turned on sometime in 
the past, and has produced the helical pattern of emission visi¬ 
ble in the left panel of Figure [2l possibly a sign of jet preces¬ 
sion / merger during that episode of activity. Using LOFAR, we 
can trace the relic plasma from that episode of activity out to 
greater distances from the core than ever before. The spectral in¬ 


dex maps of the source indicate that this phase of activity has 
resulted in an FRI-type radio galaxy (spectrum steepening from 
the core toward the edges). Then, the AGN activity has ceased 
(perhaps gradually), and after a period of time during which the 
host galaxy moved some distance to the position we observe it 
at now, it has started again, producing the less steep spectrum 
inner lobes. These lobes have clear double morphology, they are 
clearly outlined and both morphologically and spectrally distinct 
from the surrounding steeper spectrum region. At the present 
epoch, we observe an aged FRI-like large scale structure with 
an embedded restarted mdio source, s imilar to the objects an¬ 
alyzed for example by iMurgia et al.l (l201ll) . iGiacintucci et al.l 
d2007h and iGe & Ow^ (11994 . 


An alternative possibility would be that the different radio 
components (on small and large scales) belong to activity con¬ 
nected to different cores of UGC 2489 that were active at dif¬ 
ferent times. For example, a core that could have been activ e in 
the past is the core labeled E in ISchneider & GunnI (Il982h . It 
is located (in projection) to the northeast of the currently active 
AGN host core, and is co-spatial wit h the ridge of radi o emission 
which continues out to larger scales. I^iu et al.l (12014l) identify it 
as the origin of (at least) the large scale extended emission (they 
do not comment on the double lobed source morphology indicat¬ 
ing restarted AGN activity to the southwest). The steeper spec¬ 
tral index of the ridge of radio emission, indicating old plasma 
(compared to the currently active double lobes just to the south¬ 
west of it) raise the possibility that this core hosted an active 
AGN sometime in the past and is responsible for the compo¬ 
nent of radio emission that is seen, in projection, associated with 
4C 35.06. In our opinion, this is less likely; the optical core in 
question is less massive (judging by its optical luminosity com¬ 
pared to other nearby cores) than the currently active one and 
represents a less likely host of a SMBH. 

We can estimate the AGN activity timescale by using galaxy 
kinematics. The cores of UGC 2489 are embedded in a common 
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stellar envelope (probably a remnant of previous mergers). We 
assume that the core hosting the AGN is moving within the en¬ 
velope and that the AGN is intermittently active. The radial ve¬ 
locity differenc e of the AGN host galaxy and the stellar envelope 
(as reported bv ISchneider & Gunnll98^ is 320 km s“^. Assum¬ 
ing that the radial velocity is equal to the velocity in the plane 
of the sky, the projection of the section of the orbit traversed 
can be estimated by measuring the current angular separation be¬ 
tween the AGN host galaxy and the "ridge" of steeper spectrum 
emission to the northeast. This separation is 13" or 11.3 kpc at 
the distance of UGC 2489. Under the assumption that the shut¬ 
down and restart were near instantaneous, we get a lower limit 
for the elapsed time between the two episodes of AGN activity 
of around 35 Myr. It matches the duration of the AGN dormant 
p hase of some of the f ading radio sources reported in the sample 
of iParma et all (l2007l) . and also is in agreement wi th the derived 
duratio n of the inactive phases for the sample of iMurgia et al.l 
(I 2 OIII) . It is larger than the off times found (which are on the 
order of s everal Myr) for the two double-double radio galaxies 
studied bv iKonar et al.l (l2013l) . 

We have been unable to unambiguously determine a spectral 
break in the integrated source spectrum. Furthermore, the im¬ 
ages we have lack the proper frequency coverage at high spatial 
resolution to allow us to determine a spectral break and arrive at 
a radiative age estimate for the outermost emission regions. 

We can estimate at what frequency would the spectral break 
be located for the aged plasma originating from t he past activ¬ 
ity epi sode (before the AGN shutdown). Following iMurgia et al.l 
(I 2 QIII) . the break frequency yb[GHz] is related to the magnetic 
field strength in the emission region B[jiG\ and the time elapsed 
4 [Myr] since the start of the activity as 


Vb 


1 

1 + z 



£ 1/2 ^ 


where ^ic = 3.25(1 -rz)^ is the inverse Compton equivalent mag¬ 
netic field. We assume (for simplicity) that the active and shut 
down timescales of the AGN are the same, 4 = 70 Myr in our 
case. Taking Region A as an example, we can estimate the value 
of its e qui-partition magnetic field according to lGovoni & Ferettil 
(I 2 OO 4 I) . For an electron to proton energy ratio of unity and a 
Lorentz factor of the lowest energy electrons of 7 = 700, the 
magnetic field estimate is 2.8yuG, with a corresponding break 
frequency of Vb = 3.3 GHz. However, the break frequency 
value depends critically on the value of the magnetic field which 
can vary substantially. For example, for the lowest values of the 
Lorentz factor, the magnetic field can be as high as lOyuG, with 
a corresponding spectral break frequency of 380 MHz. Such a 
value for the break frequency is plausible if we look at the spec¬ 
tral shape of Region A shown in Figure [6l left panel. 

Region E has a steeper spectrum than region A. This obser¬ 
vation, (as well as the source morphology) is supported by the 
preliminary measurement of its flux density at 150 MHz using 
LOFAR (Marisa Brienza, private communication). It may mean 
that region E is older; it may actually be the oldest source region 
tha t we observe. _ 

ISchneider & GunnI (Il982h estimate that the UGC 2489 
galaxies will merge owing to dynamical friction (as shown by 
the stellar envelope surrounding the galaxies, created by strip¬ 
ping) in at most 2 Gyr, leaving us the possibility for the AGN to 
experience a series of outbursts before the final merger. 

An interesting similar example is presented by the radio 
source 3C 338, hosted by the cD galaxy (NGC 6166) in the core 


of the galaxy cluster Abell 2199. This is a steep spectrum rad io 
source (tr < -1.5). The radio morphology (iGe & Owenll 19941) is 
suggestive of a restarted source, with at least two phases of AGN 
activity, where the oldest radio plasma appears to be offset from 
an active core. In these respects, this radio source is very similar 
to 4C 35.06. However, 3C 338 is different from 4C 35.06 in that 
the older plasma in its central regions has no optical core coun¬ 
terpart, which gives a stronger support to the argument that in the 
case of 3C 3 38 the currently active core was the so urce for the 
older plasma. [Burns, Schwendeman & Whit3 (Il983l) give a thor¬ 
ough overview of radio sources hosted by cluster core galaxies 
and suggest how the motion of the AGN host galaxy can give 
rise to a radio morphology similar to what we observe. We note 
that merging of galaxy cores on timescales of 10^ yrs. can trig¬ 
ger an activity episode, which is in agreement with the derived 
dormant time for 4C 35.06. 

In these types of restarted radio sources which are located 
in galaxy clusters, there is a possibility that the ICM medium 
is dense enough near the cluster core, increasing confinement 
pressure. It can have an impact on any older outburst of AGN 
plasma, keeping it confined for longer periods of time, and pro¬ 
longing the time in which it can be detected by low-frequency 
radio observations such as we are presenting in this work. 


4.1. Presence of cold gas 


The detection of neutral hydrogen adds a new interesting com¬ 
ponent to the study of this object. 


HI in absorption has been detec ted before in clusters (for ex¬ 
ample Abell 2597 (P K S 2322 -123) 10’ Pea. Baum & Gallim^ 


1994: iTavloreta] 

[[ T^. 

Abell 

1795 (4C 26.42, 

Van Bemmel et al.l 

2012) in 

addition 

to other famous ob- 


jects in cen ters of clusters, like Hy dra A (iTavlorl [T996h and 
Cygnus A (IStruve & Conwavl l201Ql) . Although a systematic 
inventory of the occurrence of HI in clusters is not yet available, 
the statistics seems to be biased toward observing and detecting 
HI in the center of cooling core clusters. Abell 407 likely does 
not belong to this group, having a core which is not relaxed, but 
nevertheless we detect H I. 

The fact that 4C 35.06 is a restarted radio source may sug¬ 
gest a connection between the presence of gas and the restarting 
of the radio-loud phase: indeed such a relation is often seen in 
non-cluster radio sources. There, restarted radio sources show an 
incidence of HI detections higher than in other extended sources 
(ISaikia. Gupta & Konarll2007[ IC handola , Saikia & Guptal bOlOl: 
iChandola. Sirothia & Saikiall2Qdl[]Morganti et al.ll2005lk 

However, the HI profile detected in 4C 35.06 ap- 

S >ears to be broader than in typical radio galaxies 
Gereb. Morganti & Qosterlo^ l2014l) . However, it resem¬ 
bles some of the profiles found in other clusters (in particular 
Abell 2597 and Abell 1795 mentioned above) that show 
relatively broad wings in the HI absorption profile, perhaps 
a signature of the more disturbed kinematics of the gas due 
to the cluster envir o nment , as seen also in molecular gas 
by iMcNamara et al.l (1201 4 derived from ALMA data of 
Abell 1835. 


In low-luminosity radio sources with low-excitation spec¬ 
tra like 4C 35.06, the possibility of the activity being sup¬ 
ported by gas cooling from the hot galactic halo v ia e.g. Bondi 
accre ti on has been sugge s ted by several studies (ICroton et alJ 
20061: IXilen et al.l l200d i Hardc astle. Evans & CrostonI 120071: 


Balmaverde. Baldi & Capettill2008l) . Taking the core luminosity 
measured by Liuzzo et al. ( 2010l) an d following the reasoning of 
iBalmaverde. Baldi & Capettil (l2008l) . connecting the radio core 
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luminosity to the radio jet power and relating that to the Bondi 
accretion power (Pb = 0.1 Mc^), we estimate that the accretion 
rate of 4C 35.06 is on the order of 0.2 Moyr“^. 

If part of the detected HI is associated with infalling gas 
(something possible given the width of the profile that also spans 
to redshifted velocities) it could provide the required fueling 
material. The gas could be chaotically losing angular momen- 
tum in the ISM of the AGN host a nd accreting, as suggested by 
iGaspari, Ruszkowski & Qhl (l2013l) . 

5. Conclusions 

Studying the field centered on the nearby Compact Steep Spec¬ 
trum radio source B2 0258-h35A (observed during the commis¬ 
sioning of LOFAR), we have made a 61 MHz map of 4C 35.06 
(B2 0258-r35B), a radio source hosted by the core of the galaxy 
cluster Abell 407. This source has a complex morphology on 
several spatial scales. We used our LOFAR image in combina¬ 
tion with archival VLA imaging as well as VLSS, WENSS and 
other survey data to argue that the most extended emission can 
be interpreted as a remnant of a past episode of AGN activity. 
The integrated spectral index is moderately steep and shows no 
signs of spectral curvature. We confirm the existence of an ex¬ 
tended emission region and derive spectral index maps of the 
inner parts of the source closer to the AGN. 

We discuss two possible scenarios regarding the nature of 
the source and conclude that the best explanation is that we ob¬ 
serve two distinct episodes of AGN activity in 4C 35.06. Fitting 
a radiative ageing model to the spectrum of the central source 
region (encompassing the currently active AGN as well as aged 
plasma) we conclude that it is at most 10 Myr old, while age esti¬ 
mates (based on galaxy kinematics and radio spectral shape) for 
a source region farther away from the AGN show it to be older 
by at least a factor of a few. Although the outer region age es¬ 
timates are less precise, they support the suggestion of multiple 
phases of activity. 

We also present a discovery of HI in absorption at the loca¬ 
tion of the (currently active) AGN host galaxy. The HI detection 
has a broad profile, suggestive of a complex gas dynamics and a 
possible outflow. 

A better handle of the timescales involved can be obtained 
by better spectral coverage at sufficient resolution to disentangle 
the various source structures. Depending on the relative contri¬ 
butions of the smaller scales, our study shows that it is possible 
to use the integrated source spectrum as a proxy to learn more 
about the spectral behavior on smaller scales and the source ac¬ 
tivity history. This is important in studies of more distant sources 
which in general would be unresolved in future surveys. It is 
usually assumed that relic sources at high redshift should have 
very steep spectra, since we observe the exponential cutoff in 
their spectra which in our frame of reference is shifted to lower 
frequencies. However, for some of these sources, restarted activ¬ 
ity may modify their (integrated) spectral shapes such that they 
would be less steep. An example of such a case is presented by 
the integrated flux density spectrum of 4C 35.06, the source we 
have studied in this work. Our observations suggest that we have 
to be careful and possibly modify our radio selection criteria of 
high-redshift sources. 
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